SIAM J. Control Optim. 

Vol. XX, No. X, pp. 1-32 



© 2011 Society for Industrial and Applied Mathematics 



CRITERIA OF STABILIZABILITY FOR SWITCHING-CONTROL 
SYSTEMS WITH SOLVABLE LINEAR APPROXIMATIONS* 

XIONGPING DAlt 

Abstract. We study the stability and stabilizability of a continuous-time switched control system 
that consists of the time-invariant n-dimensional subsystems 

x = A 1 x + B i (x)u (3)el™,(el + and«el m >), where i e {1, N} 

, and a switching signal <r(-) : R+ — > {1, . . . , N} which orchestrates switching between these subsystems 

above, where A; e M nxn ,n > 1,N > 2, mi > 1, and where B t {-): R n -4 R" xm i satisfies the 
1 condition ||Sj(a;)|| < (3\\x\\ Vx € M n . We show that, if {Ai, . . . , Ajv} generates a solvable Lie algebra 

' over the field C of complex numbers and there exists an element A in the convex hull co{Ai, . . . , Ajv} 

in K nxn such that the affine system x = Ax is exponentially stable, then there is a constant S > 
for which one can design "sufficiently many" piecewise-constant switching signals a(t) so that the 
switching-control systems 

x{t) = A CT(t) a;(t) + B a(t) (x(t))u(t), x(0) £ M n and t £ R + 

are globally exponentially stable, for any measurable external inputs u(t) £ IR^M*) with ||«(t)|| < S. 
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1. Introduction. Let K™ be the real n-dimensional Euclidean space with an 
. inner product (•, •) which gives rise to a vector norm || • || on it. In this paper, we 

will focus on the stability and stabilizability issues for the continuous-time switched 
0^ ■ control system 

c\ 

(1.1) x{t) = A a{ t)x{t) +B CT(t) (ar(i))w (x(0) e K" and t > 0), 
, the subsystems of which are time-invariant continuous-time control systems 

(1.2) x{t) = Aix{t) + Bi{x{t))u (x(0) e W\t > and u E M m *) 

for i G {1,...,N}, where A, G K' iX ' 1 for all indices i, and a: (0,+co) {1, . . . , N} is 
piecewise constant and left-continuous having at most finite number of discontinuities 
on any finite interval of M + := (0, oo). Here n > 1,N > 2 and m, > 1 all are integers. 
We assume that the matrix- valued functions Bi(x) G JJ nxmi are continuous with 
respect to x G M" satisfying the linear growth condition: 

(1.3) \\Bi{x)\\ < (3\\x\\ Vieff" 

for each i G {1, . . . , N}, for some constant (3 > 0. However Bi{x) does not need to be 
Lipschitz continuous, not even locally, with respect to x G R". 
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For any given a(t) and u(t), the switched control system (1.1) is said to be (glob- 
ally) exponentially stable, provided that for any given initial state x £ K™, its so- 
lutions x(t) = <f>(t,xo,<j,u) with a;(0) = Xo, that are absolutely continuous in t but 
not necessarily unique because of the lack of the Lipschitz condition of the nonlinear 
terms Bi(x), are such that 

X + (x ,a,u) := Hmsup-log||</>(i,£o,0-,u)|| < 
t-y+ao t 

if x(t) is forwardly complete, i.e., x(t) may be extended on R + . 

The stability issues of such switched systems include several interesting phenom- 
ena. For example, even when all the subsystems (1.2) are exponentially stable, (1.1) 
may have divergent trajectories for certain switching signals a(t); see, e.g. [4, 27]. An- 
other noticeable fact is that one may carefully switch between unstable subsystems 
to make (1.1) exponentially stable; see, e.g. [39, 13]. As these examples suggest, the 
stability of switched systems depends not only upon the dynamics of each subsystems 
but also upon the properties of the switching signals. Therefore, the stability study 
of switched systems might be roughly divided into two kinds of problems [25]: 

(Ql) one is the stability analysis of switched systems under given sets of admissi- 
ble switching signals (all switching signals or switching signals obeying some 
constraints) ; 

(Q2) the other is the synthesis of stabilizing switching signals for a given collection 
of dynamical/control systems. 
In the present paper, the question that we are concerned with is a complex of the 
above two kinds of problems. In our context, all subsystems (1.2) are not necessarily 
exponentially stable themselves, but there exists an exponentially stable system in 
their convex hull co{Ai, . .. ,An} in W ixn . For this, we want to seek some kind of 
condition that may guarantee the existence of switching systems (1.1) that are globally 
exponentially stable; and further to describe such stable switching signals. 

To describe the switching signals that we are of interest to goal here, we need to 
introduce the classical symbolic space. Let 

(1.4) S+ = {i = (<,*)+- | i k e {1, . . . ,N} Vfc > 1} 

be the one-sided symbolic sequence space, which is compact and metrizable, endowed 
with the standard product topology. Then, there gives rise to the canonical symbolic 
dynamical system — the one-sided Markovian (forward) shift transformation: 

(i-5) e-.s+^s+- ( tfc )+~ h+ ( tfc+1 )+~. 

Observe that the shift 6 is continuous and surjective, not generally 1-to-l. For any 
vector a = (ai, . . . , ajv) € with < on < 1 and a\ + - ■ -+a.N = 1, we can naturally 
define a probability measure/distribution, written as P<j, on in this way: for any 
cylinder sets of length k > 1 

(1.6) [ii,...,i k ] := {(,€ ££\i>i = h,...,Lk = ik) , 
we have 

(1.7) Pa([ii,...,ifc]) = a il ■■■a ik , 

for every words (ii,...,i k ) G {l,...,N} fe . Then, (E+,0,V S ) is an ergodic dynamical 
system, see e.g. [30, 40]; that is to say, firstly 1?&(6~ 1 B) = P^(S) for any Borel subsets 
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B of and secondly P <;(-£?) = or 1 if the probability of the symmetric difference 
W's(BA9~ 1 B) = 0. In addition, it is easy to see that every l G is a density point 
of Pg; that is, every neighborhood of i has positive measure of Pg. 

Now, we make a convention for our convenience: To any i G ZJjy , there corresponds 
a continuous-time piecewise constant left-continuous switching signal 

(1.8) a L {-): R+ -> {1, ...,JV}; o>.(t) = t* whenever A; - 1< t < A; Vfc G N, 

where and in the sequel N = {1, 2, . . . }. Here we do not care the value a L (0), since 
the initial value x(0) = xq is given previously. 

Then, associated to the collection of subsystems given as in (1.2), there generates 
a switching-control dynamical system 

(1.1) 4 x = A a ^ t )X + B a ^ t ) (x)u, x(0) eK",teK + , and«eR m "W, 

for any i G . Our problem concerned here can now be stated as follows: 

Question 1.1. Let a = (ai, . . . , ajv) G K w be a positive probability vector; that 
is, < at < 1 /or 1 < k < N and oi\ + • • • + ajv = 1. 

(1) If the linear affine equation 

x(t) = (aiAi H h aAr^4Ar)x(t), a;(0) G R' 1 and t G K+ 

zs exponentially stable then, can one design switching signals o~ L (t) with low 
switching frequency so that the corresponding systems (1. l) t steered by cr L (t), 
are exponentially stable? 

(2) What condition can guarantee that are globally exponentially stable for 
IP '3- almost sure t G ? 

Question 1.1.(1), corresponding to the above (Q2), is to find switching signals cr L (t) to 
steer the switched systems (l.l) t globally exponentially stable. And Question 1.1.(2), 
corresponding to the above (Ql), is also one of the fundamental problems for the 
stability analysis of switched systems. Here we will give a unified positive solution to 
this question under an additional algebraic condition — solvability. 

Let A = {^4i, . . . , An} C ]R™ xn be arbitrarily given. Then under the Lie bracket 
[Ai, Aj] = AiAj — AjAi for all A i7 Aj G M nx ™, A generates a Lie algebra, write ^4la, 
over the field (D of complex numbers; that is the smallest Lie algebra containing A 
over the field <D. Letting 



a la - ■Ala, -4 L a - 



1(0) ,(0) 

^LA'-^LA 



A (l) 

> "^LA 



At-l) At-l) 
-^LA ' -^LA 



A is called solvable over <D, provided that A^ A = {0„ xn } for some integer t > 1, 
where 0„ X n denotes the zero matrix in (D nx ™. For example, abelian or nilpotent A 
implies solvable; see, e.g., [19]. 

It is well known that the Lie algebra >4la plays a very important role in the 
theory of reachability and controllability, for example, see [35, 5]. When each of the 
subsystems Ai is exponentially stable, some stability criteria, for arbitrary piecewise 
constant switching signals, of (1.1) have been developed under the solvability condition 
and an additional higher regularity of Bi(x), such as smoothness or analyticity, with 
respect to the state- variable x G K n ; for example, see [29, 17, 34, 24, 2, 28, 18]. 

Under this algebraic solvability condition, our main result obtained in this paper 
can be formulated as follows: 

Theorem 1.2. Assume that a — (0.1, . . . , ajy) G is a positive probability 
vector. Let A = {A\, . . . , An} C IR nxri be solvable over the complex-number field (D. 
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If it holds that 

(*) x(t) = (aiA\ + - ■ ■ + ajqAjq)x{t) 1 x(0) G E™ and t G R + , is exponentially stable, 

then there hold the following two statements. 

(1) For Pg-a.s. i G ZJjy , the linear switched systems 

(l.l) t ,o x(t) = A^(t)x(t), x(0) G E" and i G M + 

are exponentially stable. 

(2) Moreover, if condition (1.3) holds, then for any sufficiently small e > 0, one can 

/md a SoreZ sw&sei Wc^ rf/i Ps(W) > 1 - e and a constant 5 > suc/i 
i/iai /or eac/i t G W , the switching- control systems 

{l.l)u,e x(t) = A^(t)x(t) + B aL(t) (x{t))u{t), x(0) G E" and i G M.+ 

are globally exponentially stable, for any measurable external input u(-) with 

u(t) G E m »<-« and ||u(i)|| < 5. 
Although the statement (1) of Theorem 1.2 is a direct consequence of the state- 
ment (2) in the case e = 0, the separated formulations are convenient for our argu- 
ments later. 

Remark 1.3. Under the solvability property of A, condition (*) is not only suffi- 
cient but also necessary for the statement (1) of Theorem 1.2; see Theorem J^.2 below 
for the full details. 

Remark 1.4. Our arguments presented here imply that replacing the solvability 
of A by a more general condition that A admits a simultaneous triangularization, the 
statements of Theorem 1.2 and Remark 1.3 still hold. 

Determining whether or not a set A of matrices admits a simultaneous triangu- 
larization is itself a long studied problem; for example, see [21, 3, 33]. 

Remark 1.5. The statements of Theorem 1.2 are given in an almost sure sense. 
However, since an arbitrary point i G E^ is a density point of Pg, we can choose 
a sequence of stable switching sequence lM' G E^ such that iM' — > t as I — > +oo. 
Secondly, it is a well known fact that the Hausdorff dimension of Pa is equal to its 
entropy up to a constant multiplicator [9j. So, if the barycenter (A\ + • • • + Ajy)/N of 
the convex hull co{A\, . . . , An} is stable, then the set 

{t G E^ I x = A a i t -\X is exponentially stable} 

has the same Hausdorff dimension as the symbolic space E^ under any standard 
metrics. Thirdly, from [8], it follows that for every point i in the basin of It's, o~ L (t) 
is a stable switching signal under the (-k)- condition. 

Remark 1.6. For any A G M TIX ", it is stable iff each of its eigenvalues has a 
negative real part. So, if there exists a stable A in the convex hull co{Ai, . . . , An} in 
M. nxn , then one always can find a positive probability vector a = (otx + ■ ■ ■ + ckjv) G E w 
such that the convex combination a\Ai + ■ ■ ■ + olnAn is stable. 

Let us see our Theorem 1.2 from the viewpoint of approximations of solutions of 
the linear affine equation 

(1.9) x(t) = (a 1 A 1 -\ h a N A N )x(t), x(0) G R n and t G E+ 

by solutions of the differential inclusion 



(1.10) 



y{t) G {A x , . . .,A N }y(t), y(0) G R n and t G E+. 



Criteria of Stability and Stabilizability 



5 



There is a stabilizing switching strategy, proposed in [41], as follows: If condition 
(*) of Theorem 1.2 holds, i.e., (1.9) is stable, then from the continuous version of 
the Filippov-Wazewski relaxation theorem (cf. [14, 20]), it easily follows that to any 
initial state y(0) = £, one can find an associated switching signal ojt(t) such that the 
solution y(i, £, <7f), with y(0) = £, of the switched system 

(1.10k y(t) = K,{t)y{t), y(0) e R™ and t e R+ 

converges to zero exponentially fast; that is, it holds that 

X + (£,o- e ,0) = limsup-log||y(t,£,o£)|| < 0; 

i->+oo t 

this phenomenon is called pointwise exponentially stabilizable in [36, 37]. Yet one 
cannot claim the global stability of the above switching system (1.10)^ steered by 
such switching signal 0£ (£)! That is to say, different initial states £ may define different 
switching signals a^(t) suggested by [41]. In [38] also see [36], based on the Baker- 
Campbell- Haudorff formula J. Tokarzewski found a periodically switched signal cr(t) 
which enables the individual switched dynamical system 

(1.10) CT x{t) = A a{t) x(t), x(0) G R" and t G K+ 

exponentially stable. However, the period of <x(t) is sufficiently small and essentially 
the above switched dynamical system (1.10) CT defined by the periodic a(t) is a small 
perturbation of the stable system (1.9) there. Moreover, it is well known that the set 
of all periodically switched signals in Si is countable, O-Hausdorff-dimensional, and 
has only Pg-measure zero. 

Clearly, under the additional solvability or simultaneous triangularization condi- 
tions, our statement (1) of Theorem 1.2 presented in this paper is much more stronger 
than those mentioned above. For the nonlinear case, the relaxation theorem requires 
an additional Lipschitz condition for the nonlinear controlled part. There is a similar 
comparison if we additionally assume the Lipschitz continuity of Bi (x) with respect 
to x G W for all indices 1 <i< N. 

This paper is organized as follows. The rest Sections 2, 3, 4, and 5 of the paper 
are all devoted to proving Theorem 1.2. In Section 2, we will provide an exponential 
stability criterion for a time-dependent continuous-time linear equation whose coeffi- 
cient matrix is upper-triangular and complex, see Theorem 2.1 below. In Section 3, we 
will introduce a continuous-time symbolic semiflow by suspension of the classical one- 
sided Markovian shift transformation (S^,8). Then, we can think of our switching 
dynamical systems as skew-product semiflows driven by the continuous-time symbolic 
semiflow. Borrowing the symbolic semiflow and Lie's theorem of triangularization, we 
can apply ergodic theory to proving Theorem 1.2. We will prove the statements (1) 
of Theorem 1.2 and Remark 1.3 in Section 4 and the statement (2) of Theorem 1.2 
in Section 5. Since in our context Bi(x) may lack the higher regularity in x, the clas- 
sical Lyapunov stability theorems cannot work here for proving the statement (2) of 
Theorem 1.2. So, we will employ in Section 5 a new tool — Liao-type exponents, first 
introduced in [10] and then perfected by the recent work [11]. In fact, there we will 
prove a more general result Proposition 5.2. 

Finally we are going to conclude this introductory section with a question for 
further study. 

Question 1.7. Let A = {Ai, . . . , An} C M. nxn and a = (ai, . . . , ajv) be a pos- 
itive probability vector. If the (*)- condition is satisfied then, do the statements of 
Theorem 1.2 still hold without the solvability condition of A? 
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2. Lyapunov exponents for linear differential equations. In this section, 
we will consider a general linear differential equation 

(2.1) ±(t) =C(t)x(t), i(0)eFandteR +I 

where C(t) = [c' iJ (i)] 1<{ <n G (D nx ™, with complex elements, is a Borel-measurable 
matrix- valued function of the time- variable t on R + . We assume that the matrix 
function C(i) is bounded on R + ; that is, 

(2.2) su Pteffi+ 1| C(i) || <oo. 

It is well known that under condition (2.2), for every initial state xq G <C n there exists 
a unique solution of (2.1), written as x{t) = $(t)xo, which is defined on R + such that 
x(0) = xq. This implies that <P(t) is the principal matrix of (2.1); that is to say, 

<£(0) = Id c » (the unit matrix) and <P{t) = C(t)$(t) for Leb-a.s. tel + . 

Here Leb denotes the usual Lebesgue measure of R+. Then, the number 

(2.3) x + :=limsup-log||$(i)|| GRU{-oo} 

t— >+oo t 

is called the (maximal) Lyapunov exponent of (2.1). Clearly, for every nonzero initial 
state xo G (D", its Lyapunov exponent 

(2.4) X + (x ) := lmsup jlog||#(t)a?o|| < X + ■ 

If X + < Oj then we call (2.1) exponentially stable. According to the classical Lyapunov 
theory, see e.g. [26], X + ( x o) can take at most n distinct values for all xo S (D n \{0}. The 
basic question is: Does there hold max{x + (x ) | x G € n \ {0}} = x + ? If C(t) = C(0) 
for all t e R_|_ then the answer is positive. In general, we will see this is still true from 
Theorem 2.1 below. 

For the system (2.1), a very interesting fact is that, generally speaking, the sta- 
bility of the time-invariant systems 

x{t) = C(T)x(t), x(0) G <D n and t G R+, 

for every T > 0, cannot imply the stability of (2.1); this point is well- illustrated by the 
classical Marcus- Yamabc example. Consider the linear periodic differential equation 



x(t) 



-2 + 2 cos 2 t l-sin2t 
-l-sin2f -2 + 2 sin 2 t 



x(t) = A(t)x(t), x(0) G M 2 and t G 



one checks that, for each T G R, A(T) admits A = —1 as an eigenvalue of algebraic 
multiplicity 2 and so x{t) = A(T)x{t) are stable for all T; however, the differential 
equation admits the exponentially unstable solution x(t) = (— e* cost, e* sinf) T , where 
T means the transpose operation of a square matrix or a column/row vector. 

Particularly, we will be interested in the upper-triangular equations. For this, our 
result is the following, which is implicitly contained in the proof of the widely known 
Perron-Lyapunov regularity theorem [26] . 

Theorem 2.1. If C(t) = [c* J '(t)] G (D" x " is upper-triangular, i.e., c IJ (t) = for 
all n > i > j > 1, for t G R+, then it holds that 

X + = max{$i | i = 1, . . . , n}, 
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where 



•&i := limsup i / 9le(c"(i)) dt Vi = 1, . . . , n. 

Thus, (2.1) is exponentially stable if and only if #j < for all i = 1, . . . , n. 

Here 5He(c) means the real part of a complex number c€ C. This result shows 
that if the "time-average" time-invariant system 



1 r 

x(t) = Cx(t), x g C" and f g K + where C = limsup — / C(i) dt, 

is exponentially stable, then so is (2.1). Our condition (*), formulated in Theorem 1.2, 
is essentially an other kind of average — "spatial average". This point will be well 
illustrated in Section 4. 

Proof. We only consider here the simple, but nontrivial, case of the order n = 2; 
the general case can be similarly proved as this. 

Now, we define the matrix- valued function #^(f) = [4>ij(t)] e (D 2x2 for t S K + 
and < A < +co as follows: 



in(i) = exp(J 



c u (r)dT 



021 (*) =0, 



and 



^12(0 = y^c 12 ( S )0 22 (s)cxp(J c u (r)drJd S , 22 (t) = exp Ijf c i2 (r)dT 

It is easily seen that for any constant A g [0, +oo], the columns of ^zi(t) form a basis 
of solutions of (2.1) in the case n = 2. So, for any Z\ > 0, we have by (2.3) 

X + = limsup i log ||^(i)||. 

Therefore, we need to prove only that \ + < max{i?i, # 2 }. We next consider the two 
solutions of (2.1) described in $A.{t) and will show that 



satisfy, respectively, that 



<f>n(t) 
021 (t) 



and 02 (t) '■= 



012 (t) 
022 (t) 



and 



X + (0i) := limsup -log ||0i (t) || = ?9i 



X + (02) := limsup -log ||0 2 (t) || = tf 2 

t-t + oo t 



for some choice of the constant A. 
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In fact, from 0, = limsup t ^ +00 j J Q lHe(c Ii (r)) dr for i = 1,2, it easily follows 

that 

$1 = limsup ■j log ||<^u(f)|| and 2 = limsup j log ||<^ 2 2 (*)||- 
Then for any constant e > 0, one can find a constant K £ > such that 



I 022 (t) || 

\4>u(t)\\ 



Now, we need only to estimate X + (</>i2) := nm sup t ^ +00 \ log ||0i2(i)ll- For that, 
observe that for each e > 0, we have 

X + (0i 2 ) < limsup 1 (log || e /o ^Md-jl + i og || f c 12 ( s )0 22 (,s) e /" cll ^ds\\ 
< limsup i {log ||0n(t)|| + log |j / tf||022(s)||/||<Ms)||ds|| 



t— >+oo 



t 



1 /■* 

< X + limsup -log || / iOTee^-^+^dsl 

t->+oo £ Ja 



Here we exploit the fact that ||c 12 (t)|| < K for some constant K by condition (2.2). 
We set A = if 2 - 01 > 0, and Z\ = +00 if 2 - 01 < 0. Then, for every e > so 
small that i? 2 — 0i + £ < if 02 — 0i < 0, we can obtain that 

X + (<M < ^1 + limsup - log ^ — 1 if 2 - 01 > 0, 

and 

1 KK e ^ 2 -^ 1+£ )* 
X+(0i 2 ) <0i+ limsup -log— ^ — if 2 - 01 <0. 

i 02 - 0i + e 

Thus, no matter what A is equal to, it always holds that 

X + {d?i2) < 0i + (02 - 0i + e) = 2 + e. 

Since e > is arbitrary, there follows that X + (4>i2) < 02- So, 

X+(<£i)=0i and X + (<£ 2 ) = i?2, 

as desired. Thus, % + < max{0i,02} and then x + = max{0i,02}. This completes the 
proof of the theorem. □ 

It is interesting to note that based on Theorem 2.1, the Marcus- Yamabe phe- 
nomenon still can occur even if in the upper-triangular case. This point may be illus- 
trated by the following linear upper-triangular equation: 



x(t) 



-TTI. 



x(t) = A(t)x(t), t e M+ and x(0) £ 



one easily checks that, for each T £ R+, A(T) admits A = — yTt < as an eigen- 
value of algebraic multiplicity 2; however, this differential equation has the Lyapunov 
exponent x + — from Theorem 2.1, not less than 0; so not exponentially stable. 
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3. Symbolic semiflow. In this section, we will embed, in the manner of sus- 
pension, the one-sided Markov shift transformation 9: —> 2$ defined as (1.5) 
in Section 1 into a continuous-time semiflow. The latter is important for us, which 
plays the role of the driving system for the switching dynamical system that we are 
interested to here. 

3.1. Suspension of (£5,6, P5). First, we can define simply a continuous-time 
semiflow on the Cartesian product topological space 1$ x M + as follows: 

tp: [0,oo) x (X$ x K+) ->2$ X R+; (i,(fc,r)) i-> (i,r + t). 

We note here that R + = (0, 00). Next, we introduce an equivalent relationship, written 
as " ~ " , for x K + as follows: 

{£ = 6 k (b) if k = t' — t > 
L = e k (0 Hk = t-t'>o. 

We then define the suspension space as the quotient space 
(3.1) S(l+) =^ + xM + /~. 

Clearly, there is no loss of generality in assuming S(l$) = {(t,t) \ 1 G 1$ and < 
t < 1} with (t, 1) and (0(<-),O) identified, namely 

(3.1) ' S(l+) = 1+ x 1/ ~ 

where I = [0,1] is the unit interval. Write the elements of S(2$) or ^-classes of 
x R + as [l,t] from here on. So, S(l$) is a compact metrizable topological space, 
which is indeed homeomorphic to the compact product space 1$ x T 1 , where T 1 is 
the unit circle. Observe that i]j can further induce a continuous-time semiflow, called 
the symbolic semiflow, on S(l$) as follows: 

(3.2) e:[0, M )xS(£+)^S(£+); (t, [l, t]) h> [t, r + 1]. 

For example, 6>(§, [t, §]) = [t, 2] = [6» 2 (t), 0] for any t G 1$. 

If we identify 1$ with the 0-section [X$,0] C S(l$), then 6* may be thought of 

#,0], noting that [l$,f + fc] = 

Given any positive probability vector a — (cti, . . . , a at) € and let Leb denote 
the standard Lebesguc measure on the unit interval I with Leb(I) = 1. Then, we 

'N , 



as the 1-time map 0(1,-) restricted to [2$ ,0], noting that [S^,t + k] = [S^,t] for 
any t G R+ and k G N since 0(1$) = 1 7 + 



naturally define the product probability on the suspension space S(2$) 



(3.3) P 5 = P 5 <£> Leb 

in this way: for any Borel sets B\ C 1$ and B2 C I, we have 

Vs([Bi,B 2 ])='P s (B 1 )-Leb(B 2 ), where [Bi,B 2 ] :={[t,r]: t G B±, r G B 2 }. 
Equivalently, we have 

¥>([i,T])dP a ([t,T])= / ( L([ t ,r])dLeb(r)ldP s (0 

s(£+) Ui J 
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for any <p € C°(S(S^), ffi) (the space of continuous functions on S(E^)), where P5 is 
defined in the same manner as in (1.7) and dLeb(r) = dr. 

Lemma 3.1. For any probability vector a = (ax, ■ • ■ , &n) > 0, the corresponding 

is invariant and ergodic for the symbolic semiflow : [0, 00) x S(S^) — > S(S^); 
that is to say, (S(U^),0,Vs) is an ergodic semiflow. 

We notice here that the statement of Lemma 3.1 still holds if a is only a probability 
vector, not necessarily positive, by ignoring the letters k whenever ctk = 0. 

Proof. We notice that the (9-invariance of V& is standard. In fact, for any ip in 
C (S(£^),M) and <t < 1, we have 



<p(0(t r ))dT s 



J p([t,T + t])dLeb(r)}dP a (0 

i+t > 

^([ t ,r])dLeb(r) \ d¥ s (t) 



1 <p([i, •]) dr + j\{[6{i), •]) drj dP a (t) 
/ ¥ j([t ) r])dTdP a (0+ / / p([0(O,r])dP a (OdT 

/ p(M)drdP a (0 + [ f ^([ t) r])dP a (0dT 
Jt Jo Js+ 

v([t,r])dP a ([t,T]); 



similarly, for any tgNwe have 



V >(e(t,[t ) T]))dP s ([t,T])= / ip([l,t])6lPs([l,t]). 

This implies that Vs. is 0-invariant. So, it remains to prove only the (9-ergodicity of 
Vs.. For that, given any < p < 1, let 



' sup{|yQ)-y(Q|: tj = , 1 < j < 1} 
and 

^(£+) = {^C°(£+l): Mp<°°}- 

Elements of J%(Z!£) are referred to as Holder continuous functions. Let C^(S(S^)) 
be the set of bounded, Borel measurable, real-valued functions on S(S^). Define 
&p(S(i:£)) as the set of the functions G € £^(S(E£)) satisfying 



17(0= f G([L,r})dLeb(r)eJf p (E+). 
Jo 



Let .A/fi nv (S (£$"), 0) and A / li nv (^jv 1 be the sets of 0-invariant Borel probability 
measures on S(Z^) and ^-invariant Borel probability measures on , respectively. 
They are both nonempty compact convex sets from ergodic theory, see [30, 40]. Then 
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from [22, Section 5] (also cf. [1, 16]), it follows that there is a 1-to-l correspondence 
between M- mv (S(Efi), 0) and M inv {S^, 9), given by 



if and only if 



Gdp = 



G([ t) T])dr}d/i(t) VG e & P (S(E+)). 



S(S+) Js+ [Jo J 

Now, if Vs would not be (9-ergodic, then [40, Theorem 6.10.(iii)] follows that 
there arc < A < 1 and fix, p2 & ■M.i nv (S(S^), 0) such that 

Va = Xpi + (1 — A)^2 and pi ^ fi- 2 - 

Let p,i <r^-> /Xi G M- mv {S^ ,6) and fa <r^-> ^2 G M.- mv {S^, 9). Then, by the definitions 
of 7r* and P<s , we have 

Mi M2, Va <^ A//i + (1 - A)^ 2 and Vs <^ P<3- 

Thus, P5 = A/ii + (1 — A)/i2 and moreover P5 is not crgodic for (9, a contradiction to 
the classical theory of Markov chains. 

This completes the proof of the lemma. □ 

For a general discrete-time ergodic scmidynamics {X, 9, ix), one can similarly con- 
sider its continuous-time crgodic suspension semiflow (S{X),&,T > fi ). 

3.2. Linear skew-product semiflow. We now see the switching dynamics 
from the point of view of skew-product semiflow. 

Let A = {Ai, . . . , An} C R nx ™ be an arbitrarily given collection of real matrices. 
Then, it gives rise to the switching dynamical systems 

(3.4) x(t) = A aL(t) x{t), x(0) G R", t G R+, and t G E$, 

where the switching signals <r L : M. + {1, . . . , N} are defined in the same way as in 
(1.8) in Section 1. For any initial state a;(0) = ieo G R™, the solution to (3.4) steered 
by l G can be expressed as follows: 



x(t,x ) 



e tA -ix if < * < 1, 

e (t-k)A Lk+1 e A Lk . . . e A L1 Xq Hk<t<k + 1 and fc G N. 



Now, on the suspension space S(S^) defined in the manner as in Section 3.1, we 
define a random matrix as follows: 

(3.5) A.: S(2J+) -> R nXn ; [l, t] i-> A CTi(t) Vi G E$ and < r < 1. 

Because V s ([E+,0}) = 0, we do not care the evaluation of A. at the section x {0}. 
Since A[ t . r+fc ] = A [efc(t ^ T ] by a L (r + k) = a g k M (r) for every [t,r] G S(X$) and k G N, 
this makes sense and satisfies 

(3.5)' Ae (ti[4iT]) = A [tiT+t] =Ar t (r+t) Vr G R+ and < t < 00. 

Then, we can obtain, from the random matrix A., a linear skew-product dynamical 
system 



(3.6) ±{t) = K 0{t ^ T]) x{t), < t < 00, x(0) G R™, and [t, r] G S(Sj 



TV ■ 
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driven by the symbolic semiflow 0: [0, oo) x S(S^) — > S(S^). 

Let {^[ t ,r] M}teK+ denote the principal matrix of (3.6); that is to say, it satisfies 
&[i,r](t) = &e(t,[i,T])$[t,T](t) a t Leb-a.s. t £ R+ with <P[ L , T ](0) = Wrk the unit matrix. 
We can check the cocycle property 

(3.7) *[*,r](*i+«a) = < Z , e(t 1 ,[ t ,r])(t 2 )%r](ii) Vii,t 2 eR+. 

In fact, we need only to notice that for each t, € for any I £ Z+, 
V Al w if < t < 1, 



*m(*) 

and for any < r < 1 

*[l,T](*) 



e tA -i if < i < 1 - r, 

e [t-(fe-r)]A tfe+1 e A, fc . . . e (i-r)A tI iffc-r<£<fc+l-T and fceN. 



This implies the cocycle property (3.7). 

From the above arguments, we can easily obtain the following result. 
Lemma 3.2. Let 

e e , A : [0,oo) x S(S+) xl"-) S(S+) x K n 

be defined in the following way: 

(i,([t,r],.T)) H- ([i,T + t],# M (t)i) V* G [0,oo) and ([i,t],x) € S(X#) x 

Then, &e,A is a continuous-time semiflow; that is to say, 

&b,a(0, ([o,t],x)) = ([l,t],x), 

&e,A(tl+t 2 ,([l,T],x)) = e S , A (t2,@0,A(tl,([L,T],x))) , 

and 



&9,A{t, (['-j T ]j %)) is jointly continuous with respect to t, [t, r] and x. 

So, &e,A ^ a continuous-time linear skew-product semiflow driven by the continuous- 
time symbolic semiflow 0. 

Proof. Based on the argument above, we need only to notice here that for all 
t e [0, oo) there 0(t, [l, t}) = [l,t + t] for any [t, r] € S(Z^), and that for each k > 1, 
ti = i[, , . . , tfe = as any two points t,, t/ arc sufficiently close in E£. □ 

This lemma is important for our discussion later, because combining with the 
above Lemma 3.1 it permits us to employ the classical ergodic theorems for the sta- 
bility of the linear switching systems (3.6). 

4. Stabilizability of linear switching systems. This section will be devoted 
to proving the statement (1) of Theorem 1.2 and Remark 1.3 stated in Section 1. Let 

A = {A U ...,A N ] c R nxn 

be an arbitrary collection of real matrices. Then, it generates the linear switching 
dynamical systems 



(4.1) 



i\t) = A a[t) x{t), x(0) € R" and t e R + , 
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where a : R + = (0, oo) — > {1, . . . , N} are switching signals. For our goal, we need only 
to consider the special systems 

(4.1),. x(t) = A ai(t) x(t), x(0) el",f£ R +! and i e £+ 

where o~ L {t) is defined in the same way as in (1.8); that is, o~ L {t) = ik if k — 1 < t < ft 
for all k G N. 

We need an important classical result of triangularization for proving the state- 
ment (1) of Theorem 1.2. 

Lemma 4.1 (Lie's theorem [19]). If A is solvable over the field (D of complex 
numbers, then one can find a nonsingular complex matrix, say T G <D nx ", such that 

Ai = TAjT -1 Vie {l,...,iV} 

all are upper-triangular. 

Note that even if the matrices Ai have real entries, those of T defined by Lemma 4.1 
may be complex [19]. 

The statements (1) of Theorem 1.2 and Remark 1.3 then follow immediately from 
the following theorem. 

Theorem 4.2. Let a. = (a\, . . . ,0!jv) e R be a positive probability vector and 
assume that A = {A±, . . . , A^} C M nx " is solvable over the field (D. Then, 

(*) x(t) = (a\Ai + • • • + a]yAN)x(t), x(0) G R™ and t e R+, is exponentially stable 

if and only if for Pg-a.s. t e X^, iftese (4.1) t are exponentially stable. 

Proof. Let a = (aj, . . .,ajy) G R w be arbitrarily given as in the assumption of the 
statement. Let T be the matrix defined as in Lemma 4.1 by the solvability property 
of A so that 



(4.2) Ai = TA.J- 1 = 



Vi G {1,...,N}. 



And let 

&e,A ■ [O.oo) x S(S+) x R" -> S(i7+) x R" 

be the associated linear skew-product semiflow defined as in Lemma 3.2. We first 
notice that for any [t,r] G S(S^), the stability of the system 

(4.3) x{t) = A e{ttM) x(t), x(0) G R' 1 and t G R+ 
is equivalent to that of the switching system 

(4.4) x(t) = A aL{T+t) x(t), x(0) G 1™ and t G R+. 

In fact, (4.3) and (4.4) are the same equation from (3.5)'. We next consider the 
(maximal) Lyapunov exponent of (4.3) given by 

(4.5) X + (M) =limsupilog||% r] (i)||. 

t-S-+oo t 



14 X. Dai 

From 

(4.6) X + ([t,r]) =limsupilog||T^ [lir] (i)T- 1 ||, 

t— >+oo t 

it follows that for every [t,r] <G S(S^), (4.3) is exponentially stable if and only if so 
is the upper-triangular system 

(4.7) x(t) = A 0{tM) x{t), x(0) £ <D" and t £ R+ 
where 

(4.8) A [t , r] = l CTt(T) V[t,r]eS(25f). 

We should notice here that, under the variable transformation z = Tx £ (D™ for 

x £ W l , 

$ [hT] (t) := T$ M (i)ir 1 : (D™ -> C" Vi e M+ 
is the principal matrix of (4.7), where ^[ t . T ](i) : W — > 1™ is naturally complexified by 

$ [ltT] (t)(x + iy)=$ [ltT] (t)x + i$ [ltT] (t)y Vx + iy£€ n , 
where i = \f—\ is the imaginary unit; thus, 

limsup - log \\$[ L , T ]{t)\\ = limsup - log ||#[ t , r ](i)||. 

t— >+oo t t— >+oo t 

Let 

(4.9) i?i(M) = limsup - / Me(a^ (T+t) )di for i = 1, . . . , n. 

T^ + oo J Jo 

Then from Theorem 2.1, it follows that 

(4.10) X + ([i, r]) = max{^([ t , r]) | i = 1, . . . , n}. 
Define the qualitative functions 

(4.11) u Jt :S(S+)^R; [ t ,T]H>5le(S« (T) ) 

for each i = 1, . . . , n. It is easy to see that each uii is bounded and Borel measurable 
such that 

(4.12a) w* e ^(S(2#)) Vi = l,...,n 

and 

(4.12b) Wi(e(*,[t,r]))=Sne(5* (T+t) ) Vi e [0, oo) and [t, r] e S(i7+), 
since 0(f, [l,t]) = [i,T + i\. Thus, by (4.9) and (4.12b) we have 

(4.13) 0i(M) = limsup^ / Wi(0(t, [*, r])) dt Vi = l,...,n 
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for any [l,t] G S(S^), 

Corresponding to the probability vector a, there are associated crgodic probability 
measures P<j for the dynamics (S^,6) and V3 for (S(S^),0) by Lemma 3.1. Thus 
by the classical Birkhoff crgodic theorem [30, 40], it follows, from (4.13) and (4.11), 
that for 'Ps-a.s. [t, r] € S(S^) and for each i = 1, . . , , n, 

i?i(M)=/ w i ([t,r])dP a = / 5te(o^ (T) )dP a ([*,r]) 
(4.14) = jf^jjf We(a£ (t) )<tt}dP a (0 

<Re(a» )dP a (0; 

here the last "=" follows from cr,.(t) = ti for all < f < 1 and any t, = (tfc)tfi G 
Therefore, from (1.7) it follows that for Ps-a.s. [l,t] G S(X^), we have 



N N N N N 

(4.15) ^([ t ,r]) ^Y.MtyVaQkJ]) =J2Y, a * a ^f} = E a ^ e ( 

3=1 fc=i j=i fc=i 3=1 



a?) 



for every i = 1, . . . , n, where [k,j] C is the cylinder set of length of 2 defined by 
the word (fc, j) in the same manner as in (1.6) in Section 1. 
Sufficiency. Let condition (*) hold. Then the equation 

i = T(aiAi + • • • + ayv^w)^^, 2 G C and i G R+ 

is exponentially stable. From 

T(aiAi + ■ ■ ■ + unAn)^ 1 = aiTAiT" 1 + • • • + ajyTAjyT^ 1 

E ^s- 1 • • ■ * 
3=1 



AT 

E «3 «3"" 
3=1 



we have from Theorem 2.1 



A' 



(4.16) 



J2 a 3^(af) < Vi = l,. 
3=1 



Thus from (4.10), (4.15) and (4.16), it follows that 

(4.17) t}) < for TVa.s. [ t , r] G S(S+). 

Because for any [t, r] G S(i7^") the exponential stability of (4.3) is equivalent to that 
of (4.1) t , from (4.17) we see that for Pg-a.s. [t, r] G S(Z^), (4.1),. are exponentially 
stable. 

Finally, let 



7r:S(^+)^r+; [t,r]^t 
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Then ■k*{'P&) = P5, i.e., P5 = Vs tt" 1 - Thus, for P^-a.s. t G Z^, (4.1) t are expo- 
nentially stable. This shows the sufficiency. 

Necessity. Let (4.1) t be exponentially stable for P,j-a.s. 1 G Then, we can 
obtain that for 'Ps-a.s. [l,t] G S(Z^), (4.3) and hence (4.7) are exponentially stable. 
So, it follows from Theorem 2.1 that for pQ-a.s. [t,r] G S(17jy), i?j([i,r]) < for 
each i = 1, . . . , n. Further, by (4.13), Lemma 3.1 and the Birkhoff ergodic theorem, it 
follows from (4.15) that 

N 

^ <x/*Ke(a") < for each i = 1, . . . , n. 
i=i 

Thus, from Theorem 2.1 once again, it follows that the equation 

(4.18) z(t) =T(ai^i + --- + a N A N )T' 1 z(t), z(0) G (D™ and t G R+ 
is exponentially stable. So, the linear affinc equation 

±(i) = ( ai Ai H h aAr^Ar)x(0, x(0) G 1™ and t G R+ 

is exponentially stable, by considering initial z(0) € T(R n ) for (4.18). This completes 
the proof of the necessity. 

This thus proves the theorem. □ 

From the proof above or directly from Theorem 4.2, we can obtain the following, 
which is an important step toward the proof of the statement (2) of Theorem 1.2. 

Proposition 4.3. Assume that a. = [a\, . . . , ajv) G M. N is a positive probability 
vector. If A = {A±,...,An} C K" xrl is solvable over the complex field (D and the 
linear affine equation 

x = (aiAi + • • • + a N A N )x, x G R" and 

is exponentially stable, then for Vs-a.s. [l,t] G S(S^), 

(4.19) x(t) = A e(tj[hT]) x(t), x(0) G R" and t G R+ 

are exponentially stable. 

Remark 4.4. Up to here, we may ask if we can press forward without letup 
until we complete the proof of the statement (2) of Theorem 1.2. Unfortunately, we 
cannot. Although we have proved until now that for Vs-a.s. [t,r] G S(Z^), (4-19) are 
exponentially stable, yet to proving the statement (2) of Theorem 1.2, we essentially 
need to show that for Vs-a.s. [t,r] G S(S^), 

(4.20) x = A ( tt[hT]) x + f(x, t), x{0) G R" and t G R+ 

are exponentially stable, where |j/(a;,i)|| < L\\x\\ for sufficiently small L > 0. When 
11/(^)^)11 < -^ll^ll 1+7 f or some 7 > 0, we could conclude our desirable result from the 
classical Lyapunov stability theorem [26]. However, the classical Perron counterexam- 
ple [32] shows that the high order ^-condition is crucial. So, to overcome the trouble 
caused by the lack of such 7, we will select carefully out "sufficiently many good" driv- 
ing points [l,t] among the Vs-o-s. [l, t], using another ergodic theorem presented in 
Section 5. 
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5. Stabilizability of quasilinear switched systems. This section will be 
devoted to proving the statement (2) of Theorem 1.2 stated in Section 1, based on 
the statement (1) of Theorem 1.2 proved in Section 4. 

Let fi(x,t) G M, n ,i = 1,...,N, be measurable functions, which are continuous 
with respect to the state- variable x G R ra , such that the following linear growth 
condition holds: 

(5.1) \\fi(x,t)\\ < L\\x\\ V.tgR" 

uniformly for t G R + = (0, +oo) and i G {1, . . . ,N}, for some constant L > 0. To 
prove the statement (2) of Theorem 1.2, it is sufficient to show the following slightly 
general stability result. 

Theorem 5.1. Assume that a — (ax, . . . , ajy) G is a positive probability 
vector. Let A = {A\, . . . , Ajv} C R nx ™ be solvable over <C, and let the linear affine 
equation 

x = (ctiAi H h a N A N )x, x G R" and t G R+ 

be exponentially stable. Then for any sufficiently small e > 0, one can find a Borel 
subset WCE$ with P 5 (W) > 1 - £ and a constant 8 > such that for every l € W , 
the switching systems 

(5.2) ±{t) = A aL(t) x{t)+ f ai(t) {x(t),t), x(0) 6M" andt e K+ 

are globally exponentially stable, whenever fi(x,t), . . . ,fN(x,t) satisfy condition (5.1) 
with L < 8. 

Indeed, for proving the statement (2) of Theorem 1.2, we need only to use Theo- 
rem 5.1 with fi(x,t) = Bi(x)u(t) for reasonable external input u(t) £ R 1 ™"^'). 

Furthermore, by (3.5)' this theorem is equivalent to the following proposition. 

Proposition 5.2'. Assume that a = (a\, . . . ,ajsr) £ M. N is a positive probability 
vector. Let A = {A\, . . . , A^} C R ,lXTl be solvable over (D and the linear equation 

(*) x = (aiAi + • • • + aArAjv)x, x(0) £ R ra and t G R+ is exponentially stable. 

Then for any sufficiently small e > 0, one can find a Borel subset Z C S(S^) with 
VsiZi) > 1 — e and a constant 8 > such that for each driving point [t, r] € Z, the 
switching system 

i = A aL (r+t)X + f„,{T+t) {x, t), x(0) G R" and t G R + 

is globally exponentially stable, if fi(x,t), . . . ,f^(x,t) satisfy (5.1) with L < 8. 

Moreover, from Proposition 4.3 we easily see that this proposition follows from 
the following more general result. 

Proposition 5.2. Assume a. = (ai, . . . , a^) G M. N is a positive probability vector, 
and let A — {A\, . . . , An} C R nx ™ be solvable over (D. If the linear switched systems 

x = A a t r+ f\X, x G R" and t G R+ 

are exponentially stable for Vs-a.s. [t,r] G S(Z^), then for any sufficiently small 
e > 0, one can find a Borel subset Z C S(S^) with Vs(Z) > 1 — e and a constant 
8 > such that for each driving point [<-, r] G Z, the switching system 

(5.3) x = A a r T+t \x + f a j T +t)( x ) t), x G R" and t G R+ 



18 



X. Dai 



is globally exponentially stable, if fi(x,t), . . . ,fw(x,t) satisfy (5.1) with L < 8. 

From here on, we let A = {A\, . . . , An} be an arbitrarily given collection of real 
n x n matrices, N > 2. The remaining part of this section will be devoted to proving 
Proposition 5.2. 

5.1. Liao-type exponents and a criterion of stability. In this subsection, by 
introducing the so-called Liao-type exponent, we will provide a criterion of asymptotic 
exponential stability for a kind of deterministic switching systems that are defined by 
switching the following infinite number of non-autonomous subsystems: 

(5.4) x = Si(t)x + Fi(x,t), (i,i)el + xl"; ieN, 

where, for each index i E N = {1, 2, . . . }, S t (t) = Is^i)] E R nx ™ is contin- 

uous upper-triangular and Fi(x,t) E R" is continuous with respect to x and Borcl- 
measurablc in t, such that 

\\Si(t)x\\ < a\\x\\ and \\F^{x, t)\\ < t(t)\\x\\ Vx E R" 

uniformly for i E N, where a,£(t) both arc independent of the indices i E N. We note 
that is is not required that Fi(x, t) be Lipschitz with respect to x E R ra . 

Given any constant T* > 0, we let <x : R + — > N be an arbitrarily given T*-switching 
signal piecewise constant with a switching-time sequence : = T < T\ < T 2 < 

■ ■ ■ with Tj. — > +oo; that is to say, 

er(i) = i k e N whenever T fe _i < t < T k and T k - T k ^ <% Mk E N. 

Then, tr(t) defines a quasi-linear switching system 

(S,F) CT x(t) =S <T(t) {t)x + F (T{t) (x, t), (t,x)ER+xK n . 

To prove Proposition 5.2 we will present a criterion for the exponential stability of 
this type of switched dynamical systems here. 

5.1.1. Liao-type exponents of (S,F) . Let {k m }^ be an arbitrarily given 
integer sequence such that 

k = and 1 < k m - k TO _i < A Vra E N, 

where A is a positive integer. According to [10, 11], associated to this sequence {k m }g°, 
the real number 

1 m ~ 1 ( /" Tk *+i ] 

xt(S<r) ■= limsup — i / S'J, t Jt) dt 

is called a Liao-type exponent of (S, F) CT . 

Clearly, to different integer sequences {k m }+°° , one may get different Liao-type 
exponents for the same switched system (S,F) ct . See [10, 11]. 

From the above definition, it is easy to see that X^(S CT ) is independent of the 
"perturbation term" F = {Fi(x, t)}^. Let us see the linear switching system 

S CT z(t) = S^ t )(t)z(t), z(0) E R" and t E R+ 
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as the "linear approximation" of (S,F) . Then from Theorem 2.1, 
X + (S CT ) = max J lim sup ± / (t) dt 1 

l<j<n I T^+oo J Jo J 

is the (maximal) Lyapunov exponent of the linear system S CT . From the definitions 
above, we have x + (S CT ) < X+(S CT ). And, in general, X + (S CT ) ^ xt( S a-)- See [11] for 
an explicit example. 

5.1.2. Criterion of exponential stability. We now can formulate the expo- 
nential stability criterion for (S, F) via the so-called Liao-type exponent as follows: 

Theorem 5.3 ([11, Theorem 3.1]). Let Si(t) £ M. nxn be upper-triangular for each 
i £ N, and assume that So- has the Liao-type exponent xt (S<r) < associated to a A- 
sequence {k m }^ . Then, there exists a constant S > such that, whenever £(t) < S 
for t sufficiently large, the switching system (S, F) is globally, asymptotically, expo- 
nentially stable. 

Note: Here the constant 8 > is independent of F = {Fi(x, i)}igN but it actually 
depends upon the constants a,A,T*, and the Liao-type exponent xt- 

In addition, this criterion overcomes the lack of the 7-condition mentioned in 
Remark 4.4 in Section 4. 

Remark 5.4. We notice here that Lie's theorem (Lemma 4-1) deduces a simple 
upper-triangularization of A via a complex coordinates transformation T. For complex 
S = {Si(£)}igN, the proof of Theorem 5.3 above presented in [11], however, is invalid. 
So, we will need to pursue another natural and real upper-triangularization for A which 
having negative Liao-type exponents for sufficiently many switching signals o~ L (t). 

5.2. Frame skew-product semiflows and upper-triangularization. Let 

A = {Ai, . . . , An} C R nx ™ be any given. To use the criterion of stability in terms of 
Liao-type exponents, we will, in this subsection, introduce a natural triangularization 
for A different from Lie's theorem. 

5.2.1. Preliminaries. For each matrix A £ A, it gives rise to the time-invariant 
continuous-time dynamical system 

(A) x(t) = Ax(t), x £ R" and t £ R. 

Equivalently, it induces a continuous-time linear flow on the state space R": 

(A)' $a:Rx!"-> R"; (t, x) ^ e tA x. 

Sometimes, we identify &A{t,x) with <&A(f)x for all (t,x) £ R x R". Importantly for 
our goal, it can further induce frame flows as follows. 

Let 61,..., b„ £ R" be arbitrary n vectors. Then b = (pi, . . . ,b n ) is called an 
orthogonal n- frame of R" , provided that b forms an orthogonal basis of the vector 
space R™, i.e., (bi, bj) = and 7^ for 1 < i ^ j < n; if, in addition, ||6j|| = 1 for 
every 1 < i < n, then b is called an orthonormal n-frame of R™. Write F n and F\ as 
the sets of all orthogonal and orthonormal n-frames b of R™ inherited topologies from 
R" x ", respectively. Clearly, f\ is a subspace of F n and moreover Fn is compact, but 
F n is not. 

For the convenience of our discussion later, we denote the classical Gram-Schmidt 
orthonormalization procedure by "Ort^"; that is to say, for each collection of inde- 
pendent vectors u = (u\,... : u n ), Ort^u) = {v\,...,v > n ) £ F n is defined in this 
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manner: 



v\ = 


Ul 

INI' 






-tl 


U2 - 






\\U2 - 


(U2 




& = 


Un ~ 


n-1 

E 

k=l 




i - 


n-1 

E 


(«n,4)4l 



fc=l 



Similarly, the classical Gram-Schmidt orthogonalization procedure is denoted by "Ort" ; 
that is to say, Ort(tt) = (vi, ■ ■ ■ , v n ) € F n is defined by 



Vl = Mi, 



V 2 = U 2 - (^2,4)4' 



fe=l 



By induction on the dimension of the state space K n , we can represent the above 
Gram-Schmidt procedure in terms of matrices as follows: Under the canonical n- 
frame/basis of K™, e = (ei, . . . , e„), where 



ei 







o 
W 



we may see u = (ui, . . . , u n ) as an n x n matrix with columns u*\, . . . ,u n and also 
Ort'' (it) as an n x n matrix with columns v\, . . . , v\\ then 

(5.5) Qrtr(it) = u^" 1 , where € M. nxn is unique and upper-triangular. 
For any n-framc b = (b±, . . . , b n ) £ F n and any t € R, let 

#a(*,6) = (*a(*,6i),-..,*a(*,^)), 

which still is an independent collection of vectors, but not necessarily belong to Fn- 
So, for any b £ F n and any t € R, we let 

6) = Ort($ A (t, 6)) e F„ and $(t, 6) = Oxt\$ A (t, b)) e F„. 

One can easily check that 

(5.6) 3a:KxF„->F„ and : M x F* F„ 
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both arc flows in the classical sense. Specially, we should observe that b) and 
3a (i, b) both are smooth with respect to the time- variable tel. 

Inspirited by Liao's original work [23] for C 1 -class vector fields on compact Rie- 
mannian manifolds, we next introduce the so-called qualitative functions associated 
to A. 

Definition 5.5. The continuous functions 

n k :F*xA^R (k = l,...,n) 

given by 

d 



fi*(M) 



||ooI*(5A(t,6))|| V(6,A) € F n x A, 
t=o 



are called the "Liao qualitative functions of A '. 
Here 

coU : b = (bi, . . . , b n ) (->• b k Vfc = 1, . . . , n. 

Note here that, since <&A{t,x) is continuously differentiable with respect to t £ R, the 
above continuous functions are well defined. Since F\ is compact, all Ofc are bounded. 
From the definition, we easily have 

(5.7) log||coljfc(£i(T,&))|| = / Q k ($Ut,b),A)dt (k = l,...,n) 



for all T > and any b £ F\- This shows that these functions Q)- are closely re- 
lated to the Lyapunov exponents. In fact, limsup^^.^ j: J Q 6), A) dt is just 
a Lyapunov exponent of A. 

Under the canonical orthonormal basis e = (ei, . . . , e n ) of R™, we view z in R™ 
as a column vector with components z 1 , . . . , z n and b G f\ as an n-by-n orthogonal 
matrix with columns colib, . . . , col„6, successively; in addition, we sometimes identify 
a linear transformation of R" into itself with an n x n matrix. 

Given any orthonormal n- frame 6g r), we define, by linear extension, the linear 
orthogonal transformation of R" into itself 

(5.8) T b : R" -> R" 
in the way: e*k cohjb for 1 < k < n, such that 

n 

T b (z) = 6z := ^Vcol fc 6 and ||z|| = ||T 6 (z)|| VzGl". 
fc=i 

Now we define a family of linear transformations 

(5.9) # A)6 (t, •) = T. 1 (# A (t, T 6 (.))) : R" -► K" Vt £ R, 

where <^a(£, •) is defined by the equation (A) as before. Then there holds the following 
commutativity: 

R" * A ' b{t ''\ w 1 



( 5 - 10 ) tJ | We 



->• 
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Equivalcntly, in terms of matrices, 

(5.10)' $ A (t,b)=$ A (t,b)$ A , b (t,e) VteK. 

That is to say, 

OTt\$ A {t,b))=Mt,bWxl{t,e)- 
We now think of ^A.bit, •) as an n-by-n nonsingular real matrix. Then, ^ A \,{i) is 

upper-triangular from (5.5) and hence \P A ,b(t) is so. Clearly, fypA,b{t, m ) makes sense 
since $ A (t, •) is a smooth linear flow, and we have from (5.10) 



(5.H) 
Put 

(5.12a) 



A,3i(t,b)( a > !i M*>0) v Me 



&A,b{t, ' 



t=0 



and 
(5.12b) 

for any b £ F^- Then 
(5.12c) ^A,b(t + *i)=-R Xj? ji (t)6 )(ti) Vt,*iGR 

for any A £ A and any b £ f\- 

Definition 5.6. For any (A, b) £ A x F\, the linear equation 



(R A ,b) 



z(t) = R Aib (t)z(t), (t,z)eRx 



is called the "Liao triangularization system of A" under the moving frames 3^ (R, b). 
We will need the following basic results. 

Lemma 5.7. Given any (A,b) £ A x F,5, i/ie linear system (R A ^b) has the fol- 
lowing three properties: 

(1) Upper-triangularity: for any tgl, R A ^b(t) is upper-triangular with diagonal 
elements 

R k A k b {t) = {l k (3 A (t, b),A) for k=l,...,n. 

(2) Geometric property: for any zq £ R™, z(t, zq) is the solution of (R A .b) satis- 
fying the initial condition z(0, zq) = zq if and only if 

$ A {t,bz ) = $ A (t,b)z{t,z ) VieM. 

Particularly, ^F A b(t, •) is the principal matrix of (R Aj b) with ^^(0) = Idi». 

(3) Boundedness: there is a constant C < oo such that 

sup] l^>)l[<C 
teR [l<t/<n J 

uniformly for b £ f\- 
Proof. This lemma comes immediately from [7, Lemma 8] and so we omit the 
details here. □ 

As $ A (t, b) is an orthonormal n-frame of W l . it follows from the statement (2) of 
Lemma 5.7 that (A) has the same stability as (R A ,b) for any b £ F^. 
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5.2.2. Real upper-triangularization of A. For any A G A, we have con- 
sidered its induced dynamics in the last Subsection 5.2.1. We will now turn to the 
switching of these dynamics. That is to say, we consider the linear dynamical system 

&b,a: [0,oo) x S(£+) x R» -> S(S+) x R n ; (t,([t,T],x)) i-> ([t,r + f],* tM .](f, a;)), 

which is the linear skew-product semiflow, driven by the symbolic semiflow 

O: [0,oo)xS(2#)->S(2#); (t, [t,r]) ^ [t,r + t], 

defined in the same manner as in Lemma 3.2. We notice here that <P[ L . T ](t) is piccewise 
smooth with respect to t G K + . 

Using the Gram-Schmidt procedures as in Subsection 5.2.1, we have got that for 
all [t,r] G S(E£) with < r < 1 and t = (tfc)kS G 



(5.13) 

is such that 
S^r] (*> &) = 

and 
(5.14) 

is such that 

<r](*.6) = 



3i liT ](t, 6) := Ort($ [l>iT] (t, 6)) € F„ Vfo G F„ and t > 



>A tl (i,6) ifO<t<l-r, 
$ AL2 (t-l + T,$ ALi (l-T,b)) if l<i + r<2, 

^ fc+ i( t_ * + r '^», ( 1 »---.S4 tl (l-T ) 6)...)) if*<t + r<fc + l; 



3* iT] (t, 6) := Ort* (%, r] (t, b)) G F^< Vb G F>< 

6) iftG(0,l-T], 
^(<-l + r,5^ (l-r,6)) if * € (l,2]-r, 



As in (5.6), one can easily observe that 

(5.15) 5: [0,oo)xS(i:+)xF ra ^S(i:+)xF„; (t, ([t.r], 6)) ^ ([t.r + t],^^ 6)) 
and 

(5.16) 3*: [0,c»)xS(27+)xFjl->S(2#)xF 1 k (t, ([t, r], 6)) ^ ([,, r +t], 5*^(4, b)) 

both are skew-product semiflows, called frame skew-product semiflows, still driven by 
the symbolic semiflow on S(S^). 

Given any [l,t] € S(Z^) and 6 G F£, similar to Subsection 5.2.1, we can define a 
family of linear isomorphisms/nonsingular matrices: 



(5.17) 



% >T] , 6 (V)=T^ (ti6) # [liT] (t)T 6 (.) 



°[t.T] 



Vi G 



Then there holds the following commutativity: 



>V,t],6 



(*,•) 



(5.18) 



*[ t ,T](t.O 



Vi G 
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That is to say, 

(5.18) ' Ort"(# [t)T] (t,6))=«P [t)r] (t,6)^ ]ik (t,c) Viel + . 

It is easily seen that S^j, T ],b(^)) as an n.-by-n nonsingular matrix-valued function of £, 

is smooth at t with t + r Z. Thus, the left-hand side derivative 4r^[<,,T],b(*i *) makes 
sense at every i > 0. In addition, the cocycle/semigroup property holds: 

(5.19) % )T]i6 (t + s) =% Hii([liT])b)) (s)%, r])6 (i) VMGIR+. 
Similar to (5.12), we put 

(5.20) %,T],6(*) = {^W],b(*)} ^],b(*) V * > 0. 

for any [t, t] S S(Z^) with < r < 1 and any b € Fjj, noting that we identify [i, r + k] 
with [6» fc (0,r] for all k G N. 

Definition 5.8. TTie linear equation 

(%,t],&) = %,T],b(*M*)> ^(0) 6 R" and i e R + 

for any ([t, r],6) £ S(-EJy) x fjj, is called the "Liao triangularization system of A" 
under the moving frames 3^ T , (R_|_ , 6) . 

We will see that fRr t T .] j) is just a switching system of the subsystems (RA k ,b) 
for k= 1,...,N. 

For any [t, r] € S(Z^), there is no loss of generality in assuming < r < 1 from 
(3.5)'. We then observe that 

(5.21) a L (T + -):R+^{l,...,N} 
is such that 

(5.21)' o- L (T + t)=L k if (k- 1) -r <t < k -r and k € N, 

where i = (fc fc )fcS € 

The following two lemmas are useful for the study of the linear system (E[ l T ] t). 
Lemma 5.9. For any ([/,, r],b) € S(S^) x f\ with < r < 1, there holds that 



K-[t,r],b(*) 



^ tl ,b(*) i/0<i<l-r, 

R A L , ff" (fc-r b)(* _ ( fc - T )) if k - t < t < k + I - t and k e N 



where RA,b(t) * s as mi Definition 5.6. 

Proof. Since ^V-h (i) is the principal matrix of the switching system 

x = A a i T+t \x, x £ K" and t G R+, 

we have from (5.21)' 

*[t,r](0=*A tI (*) forO<i<l- T . 

So, from (5.14), (5.18), (5.20) and the statement (2) of Lemma 5.7, we can obtain the 
first equality of the statement. Similarly, we can prove the second equality. This thus 
proves Lemma 5.9. □ 

Lemma 5.10. For any ([t,r],b) <S S(Z^) x f\, there holds that 

R-[t,T],b(tl = K-ff'i(ti,([i.,r],6))(*2) Vii,i 2 e R+, 
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where (•, •)) * s ^ e frame skew-product semiflow as in (5.16). 

Proof. Without loss of generality, let < r < 1. We can find some k £ N such 
that (k — 1 ) — t < ti + 1-2 < k — t. In the case k = 1, from Lemma 5.9 and (5.12c), it 
follows that 

H[<.,r],b(*l + £2) = RA L1 ,b{t\ + h) = ^A L1 ^\ (ti,b)(* 2 )' 

on the other hand, for < r + t\ + t < 1 we have 
and so by (5.20) we have 



%1(*i>(M,b))(*2) 



7/7 



t=u 



then, Lemma 5.7.(2) follows that R[ tiT ] ib (ii + £2) = ^H(t I ,([i.,Tl,6))(*2)- F° r the case 
k > 2, the statement can be similarly proved by using the semigroup property. 
This thus completes the proof of Lemma 5.10. □ 

Recall that Ar l)T i is the random matrix defined by (3.5). From Lemmas 5.7, 5.10 
and 5.9, we can obtain the following important results. 

Lemma 5.11. Given any ([t,r],b) £ S(S^) x f\, the linear system (K,r t)T i j6 ) 
possesses the following properties: 

(1) Upper-triangularity: for any t > 0, the matrix R[ t . T ] ib (i) * s upper-triangular 
with diagonal elements 

B# r]|6 (t) = n*(8jU](*»&), for k = l,...,n. 

(2) Geometric property: for any zq £ R" , z(t, Zo) ^ fie solution of (R[ L , T ].b) 
satisfying the initial condition z(0, zq) = zq if and only if 

$ [hT] (t t bz )=tf ltT] (t,b)z(t,z ) Vi e R+. 

Particularly, S r [ t . T ] i b(i, •) is the principal matrix of (R[ tr ] jb ). 

(3) Boundedness: there is a constant C < 00 such that 



l<i,j<n 



t>0 



uniformly for b £ f\ and [i,t] £ S(S^). 

Proof. The upper-triangularity follows from the statement (1) of Lemma 5.7, 
Lemma 5.9 and Lemma 5.10. The geometric property follows from (5.18) and (5.20). 
Finally, the boundedness comes from Lemmas 5.10, 5.9 and 5.7.(3). This thus com- 
pletes the proof of Lemma 5.11. □ 

From the geometric property in the above statement, we can clearly see that for 
any driving point [i,r] £ S(S^), 

(5.22) x(t) = A aL{T+t) x(t), x(0) £ 1™ and t £ E+, 

is exponentially stable if and only if so is the Liao upper-triangular system 

(R[ t ,r],b) z(t) = R [hT]tb (t)z(t), z(0) £ R" and t £ R+ 
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for any/some n- frame b G F^. 

Comparing with Lie's theorem, our upper-triangular R[ L . T ],b(t) 1S real; yet the cost 
that we pay is that A becomes time-dependent. Importantly, from Lemma 5.10 we 
can define a natural skew-product system 

(5.23) z(t)=R& mitTm (0)z(t), zGR"andtGR+, ([t, r], 6) e S(S+) x rj}, 

driven by the frame skew-product semiflow ffi defined in (5.16). 

From Lemma 5.9, we can obtain the following lemma, which explains (H^i-ib) in 
terms of switching time-dependent subsystems. 

Lemma 5.12. Let ([t,r],6) G S(S^) x F^ be any given where < r < 1 and 
write bit) = 3? At, b) for all t G R+. Pui 




RA L1 ,b{t) fork = l, 

RA Lk ,b{k-i~T){t) for k = 2,3, ... ; 



and 

CT [liT ] : R + — » N; <t [( ., t] (£) = fc i/ (fc - 1) - r < i < fc - r VA: G N. 
Then, (R[ L . T ],b) * s ^ e switching system 

(5.24) z(i) = ^ ^(tXi), z(0) G R" and t G R+ 

and f[t,T] * s a T*- switching signal with T* = 1. 

We will see that (5.3) is equivalent to a perturbation of (5.24) and for "sufficiently 
many" driving points [t, t], (5.24) has negative Liao-type exponents. This will allow 
us apply Theorem 5.3 to proving Proposition 5.2. 

5.3. A spectral theorem. Let a = (ai, . . . , ajv) G M. N be a positive probability 
vector as in the statement of Proposition 5.2. Then, (S(S^), 0,V&) is an ergodic 
semiflow from Lemma 3.1. Since S(ZJ^) x f\ is compact, from the lifting lemma of 
ergodic measures proved in [7], we can easily obtain the following ergodicity result. 

Lemma 5.13. Let a = (a±, . . . , ajv) be a positive probability vector. Then there 
exists an ergodic probability measure tys on S(S^) x F\ for the frame skew-product 
semiflow 

& : [0, oo) x S(S+) x F\ -> S(S+) x f\ 

such that Vs is its marginal measure; that is to say, for any Borel set B C S(£$), 
V s (B)=y s (BxF* l ). 
Proof. Let 

tt: S(S+) x F\ -> S(S+); ([t,r},b) ^ [i, T ], 

be the natural bundle projection. It is continuous and surjective under the natural 
topologies. From definitions, there holds the following commutativity: 

S(2#) x F\ S(S+) x Fl 

J L v< g r+. 

S(S+) -^h S(2#) 
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As S(27jy ) x f\ is compact and V3 is an ergodic probability measure of on S(S^) 
by Lemma 3.1, it follows from [7] that there exists at least one ergodic probability 
measure, say tys for the semiflow & on S(S^) x F%, such that V& = CPs o 7r _1 , as 
desired. □ 

We notice here that [7] is for flows but same arguments still work for scmiflows 
and discrete-time continuous transformations. 

Based on Lemma 5.13, we can obtain the following result, which is an important 
step towards proving Proposition 5.2. 

Lemma 5.14. Under the same context as Proposition 5.2, the Liao upper-triangular 
systems (Inj^) are exponentially stable for ty^-a.s. ([/., r],b) £ S(S^) x f\. 

Proof. By the assumption of Proposition 5.2, we see that for 'P ( j-a.s. [t, r] € S(S^), 

x(t) = A au{j+t) x{t), x(0) G K" and i G K + 

are exponentially stable. Then, the statement follows immediately from Lemma 5.13 
and Lemma 5.11.(2). □ 

Noting that {\P[ L ^ b (t)}t£R + is the principal matrix of (E.[i, jT ] ) b) and tys is ergodic, 
from Lemma 5.14 we easily obtain the following. 

Corollary 5.15. Under the same context as Proposition 5.2, there exists a con- 
stant ~y~t < such that 

t J™ £ 1 °g|| S W] 1 6(*)|| =*3 

for¥ s -a.s. ([i,T],b)eS(E+)xFl 

Proof. Noting (5.23), this is a simple direct result of the classical multiplicative 
ergodic theorem [15, 31] and Lemma 5.14. □ 

Recall that from Lemma 5.11, the non-autonomous coefficient matrix of (K.[j., T ],&) 
is R[ L , T ].b{t) that is real upper-triangular with diagonal elements 

Rf* T] , 6 (t) = nfc(3f tlT] (*, b), A [4 , T+t] ) for k = 1, . . . ,n, 

where A[ t T+t ] = A CTi _( T+t ) as in (3.5). For our convenience, we introduce the following 
concept. 

Definition 5.16. For A = {A x , A N } C R nxn , the functions 

(5.25) ^.:S(£+)xr|4K; ([t, r], 6) ^ O fc (6, A„ l{t) ) 

for k = 1, ...,n, where flk * s ffl s *^ Definition 5.5, are called the "Liao qualitative 
functions" of the skew-product system (5.23). 

Then, from Lemma 5.11 there follows at once that the following holds. 

Lemma 5.17. The functions uik, k = 1, ...,n, all are bounded B or el- measurable 
on S(JCjy ) x Fn> such that 

Bf* T]ib (i) = ^^,(M,b))) 

for all t£l + and any Qi, r], b) G S(Z^T) x Fj*. 

Proof. We need only to prove the Borel measurability. Noting <j u (t) = i\ for all 
< r < 1, this measurability follows from (3.1)'. Note here that the section x {0} 
has null measure in S(E^). □ 
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Now, we can use the Birkhoff crgodic theorem again to obtain the following spec- 
tral theorem which presents an integral expression of the Lyapunov exponent xtt ■ 
Lemma 5.18. Under the same context as Proposition 5.2, it holds that 



X$ = max S / u k dq3 a | k = 1, . . . , n > . 

{JS(£+)XFH J 

Here x% * s given by Corollary 5.15. 

Proof. From the Birkhoff ergodic theorem [30, 40], we see that for k = 1, . . . , n, 



o> fe d^5= lim 1/ w fc (ffH(i,([t,r],6)))dt 



/ 



for Cp^-a.s. ([i, t],6) € S(Z^) x f\- Then, the statement follows immediately from 
Corollary 5.15, Lemma 5.17 and Theorem 2.1. □ 

Next, we will show that for <P,;-a.s. ([t, t],6) € S(Z^) x F„, (H[ t ,r],b) have got 
negative Liao-type exponents xt ■ To this end, we need an other ergodic theorem. 

Theorem 5.19 ([12]). Let <j>: [0,oo) x X — > X be a semiflow on a compact 
metrizable space X , which preserves a probability measure fi, and assume {t{}^2. 1 is 
an arbitrarily given real sequence with property 

h > 1, U +l = 2U \Ji > 1. 

Then, for any n real-valued functions fk(') € ^(I,/!),/; = 1, ...,n, there exists a 
Borel subset of /i-measure 1, write as T, such that for all x £ T, 

1 ' T 



and 



ft(x)= lim - f k m,x))At (fc = l,...,n) 

T->+oo 1 







lim { lim - V max \ \f*(x) - - / f k (cf>(t, x)) <it\\ ) = 0. 

i— too I f-i-oo £ * — * l<fc<n y ti J j t . 

Particularly, if fk(') S -^E°(^'?/ i ) arl ^ A* * s 4>-ergodic, then T is (^-invariant. 

This is a strengthened version of the classical Birkhoff crgodic theorem. We will 
apply it to the case where X = S(S^) x F^, <j> = ^r an d /&(•) = <*>&(•) for k = 1, . . . , n. 
Let 



(5.26) xt = / Wfcd^a Vfc = l,...,n. 

Js(e£)xf& 

Then, %i = max{x^ | fc = 1, . . . , n} from Lemma 5.18. And as a result of Theo- 
rem 5.19, we have the following corollary by choosing ti = 2 l ~ 1 for all i = 1,2,... 
and letting fi = Vps: 

Corollary 5.20. Under the same context as Proposition 5.2, there exists an 3^- 
invariant Borel subset r of S(S^) x f\ with *}3 ^-measure 1, such that 

{ i t± { i rU+i)^- 1 ) ) 

Km { , lim 7 E ,55? 1 I** " oi=I / (k r], 6))) dt| U = 

i— s-oo I >oo I * — J l<fc<n ^ Z J 
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for any ([t,r],6) £ -T. 

From this corollary, we now can choose the important Liao-type exponents as 
follows: 

Lemma 5.21. Under the same context as Proposition 5.2, for any e > suffi- 
ciently small, one can find an integer I > 1 and a Borel subset Z E ofS(S^) x F\ such 
that 

(1) ¥a(Z £ ) > 1 - e; 

(2) for any ([t,r\,b) £ Z e , 

1 ^ f r(i+i)2' ^1 

limsup — - V max { / « fc ($*(*, ([<-, r], 6))) dt\ <X%+£- 

Note: Here I can be sufficiently large and \ t is defined as in Corollary 5.15. 

Proof. This can be proved by an argument similar to that of [6, Theorem 3.7]. So, 
we omit the details here. □ 

Because the qualitative functions u)k are bounded by Lemma 5.17, we can improve 
the statement of Lemma 5.21 by choosing sufficiently large £, as follows: 

Corollary 5.22. Under the same context as Proposition 5.2, for any e > so 
small that Xs + £ < 0, one can find an integer £ > 1 and a Borel subset Z e of the 
driving space S(U^) x f\ such that 

(1) ¥s(Zs) > 1 - e; 

(2) for any ([t, r],6) £ Z E with < r < 1, letting {Tfe}^^ be defined by 

T Q = 0, T fe = fc-rVfc = l,2,..., 

we have 

Xt([i, t], b) := limsup — - V max \ / ^ w fe (^(i, ([t,r], 6))) dt\ 

ro-H-oo l m 2 e ~^ l <3< n [JT i2 i J 

< 0. 

So, for any ([/,, r],6) £ Z £ with < r < 1, X*([ L i T ]ib) lS just the Liao-type 
exponent of (R[ L . T ].b) associated to the switching signal er defined as in Lemma 5.12 
with the T^-switching-time sequence {T^jt^ and {k m }^ , where k m = m2 e for all 
m £ Z + and T* = 1. 

5.4. Proof of Proposition 5.2. Now we are ready to prove Proposition 5.2, 
which implies Theorem 5.1 and further Theorem 1.2. 

Proof. For any e > sufficiently small, let £ > 1 and Z e C S(i7jy) x F^j be given 
by Corollary 5.22. Given any ([t,r],6) £ Z e with < r < 1, wc next consider the 
stability of the switching system 

(5.27) x{t)=A C7b{T+t) x{t) + f^ {T+t) {x{t),t), x(0)£R n andt£M+, 

if fi(x, t), . . . , /jv(a;, t) satisfy condition (5.1) with L sufficiently small. Here, as before, 
R+ = (0,+oo). 

Let b(t) = $\^ T] (t,b) as in (5.16) and G(i) = T"^ for all t £ R+, where T 6(t) 
is defined in the same way as in (5.8). Then, G(t) is a family of orthogonal trans- 
formations and is piecewise smooth in t with 67(0) = [colib, . . . , col n 6] _1 . From the 
equation 

(5.28) x(t) = A aiiT+t) x(t), x(0) £ R" and t £ R+ 
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via the nonautonomous coordinates transformations 



(5.29) 



z(t) = G(t)x(t) W £ M 



+ 




So from Lemma 5.11, it follows that 




Write 



(5.31) 



F ai(x+t) (z,t) = G(t) ( T+t ) (G~ 1 (t)z,t). 



From (5.1), it follows that 



(5.31)' 



\\F aL{r+t) (z,t)\\<L\\z\\. 



Then from (5.27), under (5.29) we have got the equation 



(5.32) 



z(t) = R [ltT] , b (t)z(t) +F <rt(r+t) (z(t),t), z(0) £ R n and t > 0. 



Moreover, from (5.21)' and Lemma 5.12, the equation (5.32) becomes the following 
switching system 



which has the Liao-type exponent X^([t,T], b) from Corollary 5.22. 

Let C > be given by Lemma 5.11.(3) and x% < by Corollary 5.15. Applying 
Theorem 5.3 with a = C,A = 2 £ ,T* = 1, and xf = xt ([ L > T L b), there follows that 
one can find some constant 



such that (5.33) is globally exponentially stable if the constant L < 8. 
This completes the proof of Proposition 5.2. □ 

Then the statements of Theorem 1.2 hold. 

Acknowledgment. The author is very grateful to the anonymous referees for 
their insightful comments on this manuscript. 



[1] L.M. Abramov, On the entropy of flows, Dokl. Akad. Nauk. SSSR, 128 (1959), pp. 873-876. 
[2] A. A. AGRACHEV and D. Liberzon, Lie-algebraic stability criteria for switched systems, SIAM 

J. Control Optim., 40 (2001), pp. 253-269. 
[3] Y. A. Al'pin and K. D. Ikramov, Reducibility theorems for pairs of matrices as rational 

criteria, Linear Algebra Appl., 313 (2000), pp. 155-161. 
[4] M. S. Branicky, Multiple Lyapunov functions and other analysis tools for switched and hybrid 

systems, IEEE Trans. Automat. Control, 43 (1998), pp. 475-482. 



(5.33) 



z(t) = S^iMM+f^ttjW*),*), z(0) £ R" and t > 0, 



6 = S(C,xtX)>0 



REFERENCES 



Criteria of Stability and Stabilizability 



31 



A. BRESSAN and B. Piccoli, Introduction to the Mathematical Theory of Control, AIMS on 
Applied Math. Vol. 2, American Institute of Mathematical Science, 2007. 

X. Dai, Exponential stability of nonautonomous linear differential equations with linear per- 
turbations by Liao methods, J. Differential Equations, 225 (2006), pp. 549-572. 

X. Dai, Integral expressions of Lyapunov exponents for autonomous ordinary differential 
systems, Sci. China Ser. A: Math., 52 (2009), pp. 195-216. 

X. Dai, Optimal state points of the subadditive ergodic theorem, Nonlinearity, 24 (2011), 
pp. 1565-1573. 

X. Dai, Y. Huang, and M. Xiao, Almost sure stability of discrete-time switched linear 

systems: A topological point of view, SIAM J. Control Optim., 47 (2008), pp. 2137-2156. 
X. Dai, Y. Huang, and M. Xiao, Criteria of stability for continuous-time switched systems 

by using Liao-type exponents, SIAM J. Control Optim., 48 (2009/10), pp. 3271-3296. 
X. Dai, Y. Huang, and M. Xiao, Stability of time-varying nonlinear switching systems under 

perturbations, Preprint, 2010, arXiv:1109.1102v2 [cs.SY] 20 Dec 2011. 
X. Dai and Z.-L. Zhou, A generalization of a theorem of Liao, Acta Math. Sin. (Engl. Ser.), 

22 (2006), pp. 207-210. 

R. A. Decarlo, M.S. Branicky, S. Pettersson, B. Lennartson, and P.J. Antsaklis, 
Perspectives and results on the stability and stabilizability of hybrid systems, in Proc. 
IEEE: Special Issue Hybrid Systems, 88 (2000), pp. 1069-1082. 

A. Fryszowski and T. RzezuCHOWSKI, Continuous version of Filippov-W azewski relaxation 

theorem, J. Differential Equations, 94 (1991), pp. 254-265. 
H. Furstenberg AND H. Kesten, Products of random matrices, Ann. Math. Statist., 31 
(1960), pp. 457-469. 

B. M. Gurevich, Construction of increasing partitions for special flows, Theory Probab. 
Appl., 10 (1965), pp. 627-645. 

L. Gurvits, Stability of discrete linear inclusion, Linear Algebra Appl., 231 (1995), pp. 47—85. 
H. Haimovich, J.H. Braslavsky, and F. E. Felicioni, Feedback stabilisation of switching 

discrete-time systems via Lie-algebraic techniques, IEEE Trans. Automat. Control, 56 

(2011), pp. 1129-1135. 

J. E. Humphreys, Introduction to Lie Algebras and Representation Theory, GTM 9, Springcr- 
Verlag, New York, 1972. 

B. INGRALLS, E. D SONTAG, AND Y. Wang, An infinite-time relaxation theorem for differential 
inclusions, Proc. Amer. Math. Soc, 131 (2003), pp. 487-499. 

T.J. Laffey, Simultaneous triangularization of matrices — low rank case and the nonderoga- 
tory case, Linear Multilinear Algebras, 6 (1978), pp. 269-305. 

S.P. Lalley, Distribution of periodic orbits of symbolic and Axiom A flows, Adv. Appl. 
Math., 8 (1987), pp. 154-193. 

S. Liao, Certain ergodic properties of a differential system on a compact differentiable mani- 
fold, Acta Sci. Natur. Univ. Pekinensis, 9 (1963), pp. 309-327. 

D. LlBERZON, J. P. Hespanha, and A. S. Morse, Stability of switched systems: a Lie-algebraic 

condition, Systems & Control Letters, 37 (1999), pp. 117-122. 

H. Lin and P.J. Antsaklis, Stability and stabilizability of switched linear systems: A survey 
of recent results, IEEE Trans. Automat. Control, 54 (2009), pp. 308-322. 

A. Lyapunov, The General Problem of the Stability of Motion, Taylor & Francis, 1992. 

M. MARGALIOT, Stability analysis of switched systems using variational principles: An intro- 
duction, Automatica, 42 (2006), pp. 2059-2077. 

M. MARGALIOT and D. LlBERZON, Lie-algebraic stability conditions for nonlinear switched 
systems and differential inclusions, Systems & Control Letters, 55 (2006), pp. 8—16. 

K. S. Narendra and J. BALAKRISHNAN, A common Lyapunov function for stable LTI systems 
with commuting A-matrices, IEEE Trans. Automat. Control, 39 (1994), pp. 2469-2471. 

V. V. Nemytskii and V. V. Stepanov, Qualitative Theory of Differential Equations, Princeton 
University Press, Princeton, New Jersey 1960. 

V. I. Oseledec, A multiplicative ergodic theorem, Lyapunov characteristic numbers for dy- 
namical systems, Trudy Mosk Mat. Obsec, 19 (1968), pp. 119-210. 

O. Perron, Die Ordnunfszahlen linearer Differentialglwichungssyteme, Math. Zs., 31 (1930), 
pp. 748-766. 

H. Radjavi and P. Rosenthal, Simultaneous Triangularization, Springer- Verlag, New York 
2000. 

H. Shim, D.J. Noh, and J.H. Seo, Common Lyapunov function for exponentially stable 
nonlinear systems, J. Korean Institute of Electrical Engineers, 11 (2001), pp. 108-111. 

E. D. SONTAG, Mathematical Control Theory: Deterministic Finite- dimensional Systems, 2nd 

edition, TAM 6. Springer- Verlag, New York Tokyo 1998. 



X. Dai 



Z. Sun, Stabilizability and insensitivity of switched linear systems, IEEE Trans. Automat. 

Control, 49 (2004), pp. 1133-1137. 
Z. Sun, Stabilization and optimization of switched linear systems, Automatica, 42 (2006), 

pp. 783-788. 

J. Tokarzewski, Stability of periodically switched linear systems and the switching frequency, 

Int. J. Systems Sci., 18 (1987), pp. 697-726. 
V. I. Utkin, Variable structure systems with sliding modes, IEEE Trans. Automat. Control, 

22 (1977), pp. 212-222. 

P. Walters, An Introduction to Ergodic Theory, GTM 79, Springer- Verlag, New York, 1982. 
M. A. Wicks, P. Peleties, and R. A. DeCarlo, Switched controller synthesis for the 
quadratic stabilization of a pair of unstable linear systems, Eur. J. Control, 4 (1998), 
pp. 140-147. 



